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POSITION  OP  MC  MURDO  SOUND,  ANTARCTICA  TRACKING  STATION 
AS  DETERMINED  BY  ARTIFICIAL  EARTH  SATELLITES 

During  the  Antarctic  summer  (19&-I965)  a  doppler  tracking 
station,  under  cognizance  of  the  United  States  Antarctic  Research 
Program,  National  Science  Foundation,  was  installed  at  McMurdo  Sound, 
Antarctica  with  equipment  furnished  by  the  U.  S.  Navy.  Initial  station 
coordinates  were  acquired  by  taking  an  average  value  from  various 
sources.  Since  there  are  no  known  triangulation  or  Hiran  ties  from 
Antarctica  to  any  existing  datum  in  the  world,  it  was  expected  that 
the  coordinates  were  probably  in  error  with  respect  to  the  APL  Reference 
System  (Heuring,  1964) .  This  was  soon  verified  after  some  passes  of 
data  were  processed  by  our  tracking  programs.  Herein  will  be  described 
the  method  and  results  of  refining  the  position  of  the  McMurdo  Sound 
Doppler  Tracking  Site  (hereinafter  referred  to  as  Station  019) . 

Method 

Prior  to  defining  the  details  of  the  calculation  a  brief 
foundation  will  be  laid.  Passes  of  doppler  data  are  received  from 
stations  in  the  TRANET  doppler  tracking  system  (for  an  overall  discussion 
of  the  tracking  system  see  Newton,  1963) .  These  data  are  processed  and 
compared  to  similar  data  generated  by  an  updated  orbit,  i.e.,  one  derived 
from  integrating  the  equations  of  motion  to  some  particular  time,  t, 

-4 

thereby  defining  the  satellite  position,  r(t),  and  velocity  r(t).  Then 
the  difference  between  the  experimental  and  generated  data  is  minimized 
by  adjusting  the  initial  conditions,  i.e.,  the  orbit  parameters.  At  the 
minimum,  the  orbit  is  deemed  an  improved  orbit. 
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These  calculations  are  described  in  the  orbital  improvement  programs 
(see  Holland  and  Yingling  ri963],  Black  [1963a],  Welty  [1963a 1, 

Welty  [1963b!,  Datatrol  Corporation  [1963],  Holland  [1963!,  Albright  T 1963"! 
and  Monahan  [1964]).  Using  the  improved  orbit,  an  assessment  of  the 
individual  passes  used  to  improve  the  orbit  is  made  on  a  routine  basis. 
Further,  passes  which  have  not  been  used  to  improve  the  orbit  may  be 
assessed..  The  latter  type  of  computer  run  is  termed  an  Editor  run  (ED) 
whereas,  the  former  is  a  tracking  run  (T).  Both  type  of  runs  were 
employed,  in  adjusting  the  position  of  Station  019.  Some  of  the 
tabulated  results  which  were  used  from  these  runs  are: 

(l)  along-track,  e  j 

ai 


(2) 


slant  range  error. 


(3)  satellite  elevation,  E^ 


(4)  satellite  azimuth  clockwise  from  north, 

All  of  these  values  are  evaluated  at  the  point  of  closest  approach  of 
the  satellite  to  the  station  where  the  subscript  i  refers  to  a  particular 
pass  in  a  run. 

Along-track  error  is  the  difference  in  position  of  the  station 

between  that  defined  by  the  given  station  coordinates  and  that  defined 

by  the  doppler  data  recorded  at  the  station.  Its  positive  direction  is 

that  of  the  relative  velocity  of  the  satellite,  r,  and  station,  r^,  i.e., 

•  • 

(r  -  ?T)  [Black,  1963b].  The  slant  range  error  is  similarly  defined 

except  that  the  positive  direction  is  from  station  to  satellite  [Black, 1963b!. 
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An  orthogonal  triad  is  defined  by  e  ,  e_  and  e  x  e  ,  but  the 

pi  ai  pi  ai 

magnitude  of  the  latter  is  not  evaluated  because  of  a  very  ill- 
conditioned  matrix  ("Black  1963b].  Further  details  of  this  and 
other  coordinate  systems  and  inter-relations  may  be  found  in  Guler 

r 19631. 

Contained  in  Figure  1  are  four  possible  geometric 
configurations  for  a  pass  of  a  satellite  where,  at  closest  approach, 
representation  (a)  is  a  satellite  moving  northward  to  the  east  of  a 
station,  (b)  is  a  satellite  moving  northward  to  the  west  of  a 
station,  (c)  is  similar  to  (a)  but  moving  southward  and  (d)  is 
similar  to  (b)  but  moving  southward.  Construct  a  right-handed 
cartesian  coordinate  system  whose  X  and  Y  axes  are  in  the  station 
plane,  with  +Z  in  the  station  vector  direction  out  of  the  earth 
and  +X  in  the  direction  of  the  projection  of  the  positive  slant 
range  vector  onto  the  station  plane  (Figure  le).  In  this  X,  Y,  Z 
coordinate  system  we  will  express  the  displacement  of  the  position 


of  the  station  from  its  initial  position  according  to  e  ,  e  ,  A. 

ai  pi  1 


and  First,  determine  X^  Y^,  Z^. 


N  N 


«o  ALONG  TRACK  ERROR 
>p  SLANT  RANGE  ERROR 
E  ELEVATION 

A  AZIMUTH  FROM  NORTH 
A  TRACKING  STATION 
O  SATELLITE 
X,Y  STATION  PLANE 
Z  STATION  ZENITH 


Z  (•) 


FOR  (b)  AND  (c),  «0  IS  SET  EQUAL  TO  -i0  BEFORE 
COMPUTATIONS  BEGIN,  SEE  TEXT. 


Fig.  1  GEOMETRIC  CONFIGURATIONS  FOR  A  SATELLITE  PASS 
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the  +  Y  direction,  whereas  in  Figure  1  (b,c)  +  e  is  in  the  -Y  direction. 
For  the  latter  cases,  i.e..  Figure  1  (b,c),  set  e„  =  -  e  so  that  all 

J  70  f  ft 

e  >0  are  in  the  +Y  direction.  Ihe  adjusted  station  position  is  X^ 

It  remains  to  express  ,Y^;Z^  in  the  APL  Reference  System,  i.e.,  X^,  YA, 
Z^.  The  latter  is  here  accomplished  by  three  rotations  and  a  translation. 
First  rotate  the  X,  Y,  Z  coordinates  about  the  Z  axis  by  \|r^  =  A.^  -  rr/2  so 
as  to  express  the  new  position  in  a  coordinate  system  X',  Y*,  Z'  with  X1 
toward  the  east  and  Y*  toward  the  north,  i.e., 


/X,\  = 

/  \ 

/  , 

1 Yi  - 

cos  \|r. 

! 

sin  \|^ 

°\ 

\ 

xi 

-  sin 

i 

cos  \|ri 
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\ 

V  o 

0 
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i  = 


where  I  is  the  number  of  passes  used  to  adjust  the  station  position. 

i  »» 

Rotate  about  the  X  axis  such  that  the  Z  axis  is  parallel  to  the  Z^  axis, 
i.e.,  C  =  -  n/2  where  <p,p  is  the  station  spherical  latitude  used  in  ED. 


X  \  =  /I 
1  \ 

\ 


Y.  '=  0 

1  i 


cos  C  sin  C 


i  =  1,  . . .,  I 


,0 


sin  C  cos  C 
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"  at  m 

Rotate  about  the  Z  axis  such  that  the  X  and  Y  axes  are  parallel  to 
and  axes,  i.e.,  B  =  -(Xq,  +  tt/2)  where  X,p  is  the  station  spherical 
longitude  used  in  ED. 
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X  «  tan”1  E/x] 
r  =  [X2  +  y2  +  z2]* 

With  these  coordinates  an  ED  is  made  to  test  the  validity  of  the  station 

adjustment,  i.e.,  determine  if  the  rms  total  error  (e  and  e  )  of  the  ED 

pi  ai 

has  been  reduced.  Ciis  process  is  continued  until  the  rms  total  error 
is  less  than  that  of  the  orbit. 

Results 

Tracked  orbits  (day  043  and  day  U5 2)  of  satellite  03164,  defined 

with  no  Station  019  data,  were  used  with  the  Editor  programs  yielding  slant 

range  and  along-track  errors  which  allowed  for  adjustment  of  the  Station 

019  coordinates.  Table  1  contains  the  step-by-step  results  where  data 

span  indicates  the  period  in  which  the  data  used  from  Station  019  were 

observed;  o^,  is  the  rms  error  of  the  tracked  orbit;  c^  is  the  rms  error 

of  the  Station  019  data  from  an  Editor  run;  Station  019  passes  refers  to 

the  number  of  passes  used  in  an  Editor  run;  ED  coordinates  are  those  used 

in  a  particular  Editor  run;  and,  the  adjusted  coordinates  are  the  result 

of  correcting  the  ED  coordinates  via  the  Editor  run  results,  i.e.,  e  e 

pi  ai 

etc.  3he  adjusted  coordinates  become  the  ED  coordinates  for  the  next 
Editor  run  if,  fo~  a  particular  run,  >  o^.  If  <  oT,  further  Editor 
runs  with  the  same  orbit  are  meaningless,  for  the  station  position  error 
is  less  than  the  orbit  error.  This  is  the  situation  after  runs  four  and 
six,  and  since  the  specified  are  typical  values  for  standard  tracking 
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runs  the  adjusted  coordinates  of  run  six  would  be  deemed  the  best 
coordinates  for  Station  019 . 

The  validity  of  these  coordinates  was  assessed  by  tracking 
satellite  63041  on  day  052  -  1965  with  nine  passes  of  Station  019 
data  and  64  passes  of  data  from  Stations  003,  008,  Oil,  012,  013,  014, 

018,  106,  111,  and  115  whose  coordinates  are  well  defined  in  the  APL 
Reference  System  (see  Ouler  [1965 ]) . 

Table  2  contains  the  tracking  results  using  the  run  six 
adjusted  coordinates  for  Station  019 .  The  a T  for  all  passes  used  in  the 
tracking  is  greater  than  the  aT  for  the  Station  019  passes;  thus,  the 
adjusted  coordinates  (run  six  Table  l)  are  considered  a  good  representation 
of  the  Station  019  position. 

Using  the  x,  y,  z  values  (Eq.  l)  for  the  adjusted  coordinates 
of  run  six  in  Table  1  an  approximation  can  be  made  for  the  geodetic 
coordinates  <p^,  Xq  in  the  APL  Reference  System  (Heurlng  [1963])  and 
appears  in  Table  3  along  with  the  initial  station  coordinates.  This 
is  an  approximate  value  because  the  distance  h  from  x,  y,  z  to  the  APL 
Ellipsoid  (a  =  6378. 166  km  and  f  =  1/298. 3)  includes  the  geoidal  height 
as  well  as  the  station  elevation  above  mean  sea  level  (ostensibly  the 
geoid)  with  only  the  latter  known.  However,  the  station  elevation 
(35.6  m.)  and  distance  h  (~  -33.5  m,  implying  the  ellipsoid  is  above 
the  point  x,  y,  z)  yield  a  geoidal  height  value,  H,  of  about  -70  m 
which  is  close  to  the  geoidal  height  value  of  -50  in  Guier  1*1965].  In 
view  of  the  small  values  h,  5  and  H  the  error  incurred  in  specifying 
9  ,  \q  from  x,  y,  z  is  of  the  order  of  a  meter  or  less.  It  is  not 


_ _ TABLE  3 

INITIAL  AND  FINAL  GEODETIC  C 
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meant  to  imply  that  the  position  is  known  to  a  meter  only  that  the 
transformation  from  spherical  to  geodetic  coordinates  is  specified 
with  this  accuracy,  The  accuracy  of  the  position  depends  on  the 
orbit  accuracy  which  in  thin  case  is  less  than  70  meters  rms*-wise. 

This  report  is  considered  preliminary  being  based  on  a 
email  data  sample  sent  to  us  by  the  last  possible  mail  prior  to 
their  1965  winter.  Accordingly,  a  more  thorough  study  is  contemplated 


at  a  later  date . 
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